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Abstract 

A one-dimensional driven diffusive system with two types of particles and nearest 
neighbors interactions has been considered on a finite lattice with open boundaries. 
The particles can enter and leave the system from both ends of the lattice and there 
is also a probability for converting the particle type at the boundaries. We will show 
that on a special manifold in the parameters space multiple shocks evolve in the 
system for both species of particles which perform continuous time random walks 
on the lattice. 
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1 Introduction 



Recently much attention has been paid to the study of one-dimensional driven 
diffusive systems because these relatively simple systems show a variety of crit- 
ical phenomena such as out-of-equilibrium phase transitions and spontaneous 
symmetry breaking [1,2]. Another remarkable feature of these systems is that 
under special circumstances, for example by imposing restrictions on micro- 
scopic reaction rates, shock waves evolve in the system in the steady state [3]. 
A simple model of this kind is the Totally Asymmetric Simple Exclusion Pro- 
cess (TASEP). The TASEP is a system of particles with hardcore exclusion 
interactions which hop only in one direction on a one-dimensional lattice with 
open boundaries. This diffusion model can describe hopping conductivity in 
ionic conductors, traffic flow and interface growth [1,2]. The TASEP shows 
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rich non-equilibrium behaviors such as shock wave [4,5,6], boundary induced 
phase transition [7], and dynamical scaling in the universality class of the 
Kardar-Parisi-Zhang equation [8,9]. This model has been solved exactly using 
a so-called Matrix Product Formalism (MPF) [10]. According to this formal- 
ism the stationary state probability distribution function of the system can 
be written in terms of expectation values of products of two non-commuting 
matrices which are associated with the existence of particles and holes at each 
site of the lattice and satisfy a quadratic algebra. 

Another model of this kind is Partially Asymmetric Simple Exclusion Pro- 
cess (PASEP) with open boundaries. In this model, which is defined on an 
integer lattice Z of length L, each site of the lattice i (1 < i < L) is either 
empty Tj = or occupied by at most one particle T{ = 1. The system evolves 
according to a stochastic dynamical rule: during each infinitesimal time step 
dt the transitions allowed for the bond (i, i + 1) are (10 — > 01) with rate 1 
and (01 — > 10) with rate x. The parameter x is positive and measures the 
strength of the driving field. Without losing generality one can assume that 
x < 1. Particles are injected from the first and the last sites of the lattice 
with rates (1 — x)a and (1 — x)S. They also leave the lattice from both the 
first and last sites with rates (1 — x)j and (1 — x)/3 respectively. It is known 
that the phase diagram of the PASEP has three different phases depending on 
a, (3, 7, 5 and x: in the steady state it has a high-density phase, a low-density 
phase, and a maximal-current phase. Recent investigations show that for the 
PASEP with open boundaries (and even on an infinite lattice) a travelling 
shock with a step-like density profile might evolve in the system provided that 
microscopic reaction rates are tuned appropriately [11,12]. The density of par- 
ticles on the left and the right hand sides of the shock position is a function 
of these rates. The shock position then performs a random walk in the bulk 
of the system and also reflects from the boundaries. There is also a possibility 
for the existence of n consecutive shocks in the system. In an infinite system 
consecutive multiple shocks evolve according to n-particle dynamics [11]. By 
slightly different definition of the shock it has been shown that the same phe- 
nomenon takes place for the PASEP on a finite lattice with open boundaries 
[12]. The PASEP has also been studied using the MPF [13,14]. It has been 
shown that its quadratic algebra has an n- dimensional matrix representation 
[15,16] provided that exactly the same constraints necessary for the existence 
of n consecutive shocks in the system hold [12]. 

Despite the remarkable work on the dynamics and structure of shocks in the 
systems with one species of particles, not much is known about the shocks in 
the systems with more than one species of particles. In [17] the authors have 
introduced a one-dimensional driven lattice gas with two types of particles 
and nearest neighbor hopping. Part of their work is devoted to the study of a 
single shock dynamics in the system. They have investigated the time evolu- 
tion equation of a product shock measure and shown that the shock is stable 
and performs a random walk provided that some constraints are satisfied. In 
this paper we study the same model; however, we investigate the possibilities 
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for the existence of multiple shocks in the system. It turns out that by impos- 
ing some conditions on the microscopic reaction rates, stable multiple shocks 
evolve in the system. For the special case where only a single shock exists, our 
results converge to the ones obtained in [17]. We will also show that under 
these constraints the quadratic algebra of this model can be mapped to the 
one which appears for the PASEP. 

This paper is organized as follows. In section 2 we will briefly review the 
PASEP from the MPF point of view. In section 3 we will define the model 
and apply the MPF to study its steady state properties. In section 4 the con- 
ditions for the existence of multiple shocks in our model will be studied. The 
concluding remarks will be given in the last section. 



2 The PASEP 



In this section we will first review the steady state properties of the PASEP 
with open boundaries defined in section 1, because as we will see the steady 
state properties of our two-species model are closely related to those of the 
PASEP. The exact steady state properties of the PASEP for all values of the 
boundary and bulk parameters might be calculated by using the MPF [14] 
which involves the representation theory of a quadratic algebra equivalent 
to a g-deformed harmonic oscillator algebra [13]. According to the MPF the 
stationary probability distribution function of the PASEP is given by 

P({ n , • • • , r L }) oc (W\ fl( Ti D + (1 - n)E)\V) (1) 

i=i 



provided that the operators D and E besides the vectors (W\ and \ V) satisfy 
the quadratic algebra [16] 



DE-xED = (1 -x)(D + E) (2) 
(f3D-5E)\V) = \V) (3) 
(W\(aE--yD) = (W\. (4) 

The non-commuting operators D and E stand for the existence of a particle 
and a hole at each site of the lattice. It has been shown that the associated 
quadratic algebra (2-4) has exactly one n-dimensional irreducible representa- 
tion for any finite n provided that the following constraint is satisfied by the 
bulk and the boundary rates [15,16] 

x 1 - n = K + (a,7)K+(P,S) (5) 
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Fig. 1. A typical Bernoulli shock measure for n— 1 consecutive shocks for the PASEP. 
The density of particles at the shock positions h L is 1. 



in which we have denned 

+ v + 1 + \J(u - v - l) 2 + Auv 



K+(U,V) 



-u 



2u 



(6) 



Typical n-dimensional representations for (2-4) are given in [15] and [16]. 
The finite dimensional Fock representations allow us to drive exact results for 
the PASEP on some special curves of the phase diagram. Theses curves pass 
through the high-density and the low-density phases but not the maximal- 
current phase. 

The evolution of product shock measures in the PASEP has been studied both 
for open boundary condition and on an infinite lattice. In the following we 
study the dynamics of travelling shocks in the PASEP with open boundaries 
briefly. We define a Bernoulli measure for n— 1 consecutive shocks as a product 
measure with density 1 at the shock positions hi and intermediate densities pi 
between sites and fcj. This product measure can be written as follows 

|&4, fc 2 , • • • , k n -i) = 

|1) (g) |p 2 )®(fe2-fcl-l) ® |1) ® . . . ® |1) (g) |p n )®( L - fe "-l). 



in which < p\ < pi < • • • < p n < 1. A typical shock measure is plotted in 
Fig. 1. Now the question is that how (7), as an initial configuration, evolves 
in time. This is given by the master equation. It is known [11,12] that for the 
PASEP with open boundaries the time evolution equation of the Bernoulli 
shock measure (7) has a closed form and is similar to the time evolution 
equation of n — 1 random walkers on a finite lattice provided that we have 



Pi 



Pn 



1 

1 + «+(a,7) 
1 + «+(/?, 5) 



(8) 
(9) 
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and that the consecutive densities are related through 



-pi) 



P<(1- 




in which i 


= !,-• 


Pn(l - 


-Pi) _ 



Pl(l -Pn) 



1 



(10) 



(11) 



-x 



The ith random walker hops to the left and the right with rates Li = —_ p 
and Ri = respectively. Using (8) and (9) one can rewrite (11) in terms of 

K + (a, 7) and k+(/3, 5) to see that it is simply the condition (5) i.e. the necessary 
condition for the existence of an n-dimensional representation for (2-4). One 
can now construct the stationary state of the PASEP in terms of a linear 
superposition of Bernoulli shock measures with n — 1 shocks; therefore, the 
n-dimensional representation of the quadratic algebra describes the stationary 
linear combination of shock measures with n — 1 consecutive shocks. 



3 Two-species model with open boundaries 



The model is defined on a finite lattice of length L with two species of particles. 
Each site of the lattice is either empty or occupied by a particle of type A or 
B. We assume that in the bulk of the lattice the particles hop to the left and 
right according to the following reaction rules 



Bulk: 



AO -> OA with rate 1 
OA -> AO with rate x 
BO -> OB with rate 1 
OB -> BO with rate x 
AB -> BA with rate y 
BA -> AS with rate y. 



(12) 
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At the boundaries the particles are injected and extracted. The particle type 
can also be converted there. The boundary processes are defined as follows 



Boundaries: < 



— > A with rate ai( T )A 
— > B with rate a^ r ) B 
A — > with rate /?z( r )A 
£> — > with rate 
A^> B with rate 7z( r )A 
B ^> A with rate 7;( r )_B 



(13) 



in which the indexes / and r indicate the left and the right boundaries respec- 
tively. In the following we will investigate the steady state properties of this 
model using the MPF. 

We assign the operators E, A and B to an empty site, an A particle and a B 
particle respectively. By applying the standard MPF [10] we find the following 
quadratic algebra for our model defined by (12) and (13) 



AE - xEA = aE + (a + b)A 
BE - xEB = bE+{a + b)B 
y(AB - BA) = aB-bA 

for the bulk operators and also 



(14) 
(15) 
(16) 



[{fir a + 1va)A - a rA E - lrB B - a]\V) = 
[(PrA + lr B )B - a rB E - lrA B - b]\V) = 
(W\[((3 lA + Iia) A - a lA E - llB B + a] = 
{W\[(/3 lB + llB )B - a lB E - llA A + b] = 



(17) 
(18) 
(19) 
(20) 



for the boundary terms in which a and b are non-zero arbitrary numbers. 
We have looked for finite-dimensional representations of the algebra. We have 
found that a one-dimensional representation for the algebra exists in which the 
operators A, B and E are replaced by non-zero real number a, b and e. This 
representation is associated with a uniform density for each type of particles 



on the lattice. Defining these densities as p A 



and p B 



the 



a+b+e I a • a+b+e 

necessary conditions for the existence of a one-dimensional representation are 



Pa = au(l - Pa - Pb) - (Pia + 1ia)pa + IibPb 

= (PrA + lrA)pA - Oi rA (l - p A - p B ) - J rB p B 

= (1 - x)(l - p A - Pb)pa 



(21) 
(22) 
(23) 
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Pb = «ib(1 - Pa- Pb) - {PlB + lw)pB + IiaPa 
= (PrB + 1tb)Pb - a rB (l - p A - Pb) - %-aPa 
= (1 - x)(l - pa - Pb)pb 



(24) 
(25) 
(26) 



in which pa '■= i , and n B := — |— . We have also found that the algebra 

a+b+e < D a+b+e & 

(14-20) has a finite-dimensional representations provided that the matrices A 
and B commute with each other but not with E. This means that in the steady 
state the probability for finding configurations consisting of blocks of mixtures 
of A and B particles surrounded by holes does not depend on the exact position 
of the particles in each block, instead only the number of them in each block 
will be important. For instance one can consider the case A = r B in which 
r is a non-zero real number. This condition indicates that the density of A 
particles on each site is always r times the density of B particles on the same 
site. By defining a new operator D := A + B = (1 + -)A associated with the 

total density of particles on each site of the lattice and choosing a := it 
can easily be verified that (14) and (15) both converge to (2) while (16) gives 
b = j^. By defining two new parameters 5 and (3 and imposing the following 
constraints on the right boundary rates 



1 (1+r) 

S:=- (a r A + otrB) = —r, r«ri (27) 

1 — x r(l — x) 

X 1 

P ■= Z (PrA + It A — ) = Z (PrB + IrB ~ HVa) (28) 

1 — x r 1 — x 

in which a r A = r a rB one can see that the equations (17) and (18) become 
identical to (3). On the other hand one can define two new parameters a and 
7 and impose some constraints on the left boundary rates 



1 (1 + r) 

ol:=- (a tA + a w ) = —- zO-ia (29) 

1 — x r(l — x) 

7 := — !— (Aa + HA ~ —) = -r^—(PiB + iw ~ r llA ) (30) 
1 — x r 1 — x 

in which ol\ A = r a\ B to see that the equations (19) and (20) become iden- 
tical to (4). Therefore, the operator D associated with the total density of 
particles of kind A and B on each site of the lattice besides the operator E 
satisfy the PASEP quadratic algebra, provided that the constraints (27-30) 
are satisfied. Now one can easily see that on an special manifold defined by 
the aforementioned constraints in the parameters space of our model an n- 
dimensional representation exists for the quadratic algebra (14-20) provided 
that (5) is also held. To conclude, we have shown that under some conditions 
the quadratic algebra of our two-species model defined by (12) and (13) has 
exactly one n-dimensional irreducible representation for any finite n. Having 
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a finite-dimensional representation for the algebra one can easily calculate the 
physical quantities such as the mean density of particles at each site. 



4 Dynamics of Multiple shocks 



Now we investigate the shock formation and also the shock dynamics in our 
two-species model. Since the quadratic algebra of the model in terms of the 
total density of particles operator D is exactly the one for the PASEP and 
therefore, has a finite-dimensional representation, one can conclude that mul- 
tiple shocks might evolve in the system provided that the constraints (8), (9) 
and (10) obtained for the PASEP also hold for the total density of particles in 
our model defined as p := pA + Pb- However, since the density of A particles 
on each site is always r times the density of B particles at the same site it 
turns out that the multiple shocks structures exist for both densities pa and 
Pb separately. 

In the following we will study the case where there is only a single shock in 
the system in details. This will be associated with the existence of a two- 
dimensional representation for the quadratic algebra. We define the total 
density of particles on the left- (right-) hand site of the shock position as 
Pi(r) '■= pi(r)A + pi(r)B- Since A = r B we always have 



PjA = P_rA =r (31) 
PlB PrB 

r 

Pl(r)A = 7^7^M ( 32 ) 

Pl(r)B = Y^. Pl ^y ( 33 ) 

From our discussion in the previous section and using (11) we see that in order 
to have a single shock the total density of particles on different sides of the 
shock should satisfy the condition 

= i (34) 
pi{l-p r ) X 



On the other hand the total density of particles on each side of the shock 
should be obtained from (8) and (9) 



Pi 



1 

1 + K+(a,7) 



(35) 
(36) 
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in which a, (3, 7 and 5 should be replaced from (27-30). It is not difficult to 
verify that using (27-30) the condition (35) can also be written as 



a u = Pu((l -x) + P \ A + 7L4 ~~ ) (37) 

1 - PlA - PlB 



olib = Pw((l -x) + 
and (36) as 



Pib + Jib - nu , 
1 - Pia - PlB 



arA = -x) + ^ A + lrA 7;B ) (38) 

J- — PrA — PrB 

( /-, \ , PrB + IrB ~ r~f rA , 
«rB = p r s(-(l - X) H J. 

As we mentioned in section 2 the shock position hops to the left and to the 
right with the following rates: 



' " (39) 



r 



R= 1 —^. (40) 
1 - Pi 

It is also known that for the PASEP with a single shock the shock position 
reflects from the left boundary with rate L — (1 — x)(y + pi) and also from 
the right boundary with rate R = (1 — x)(^ — p r ). For our two-species model 
they become 



r a r A + a rB I / 1 \ , A ~s 
L = + (1 - x)pi (41) 

Pr 

R = (1 — x)p r . (42) 

Pi 

In [17] the authors have introduced a two-species model with open boundaries 
similar to the one studied here and investigated the time evolution of a single 
shock in the system by introducing a product shock measure. They have found 
that under exactly the same conditions which we introduced in (27)-(38) a 
stable single shock evolves in the system which performs a random walk on 
the lattice. However, in this paper we have proven that not only one, but 
also multiple stable shocks might evolve in the system under slightly different 
conditions. The expressions which we have obtained for the hopping rate of 
the shock positions in the bulk of the lattice and also the reflection rates 
from the boundaries are quit in agreement with those obtained in [17]. One 
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can easily check that the necessary conditions for having a stationary product 
measure with uniform densities for both kinds of particles obtained in [17] are 
exactly the conditions for the existence of a one-dimensional representation of 
the quadratic algebra (14-20) obtained in (21-26). 



5 Concluding remarks 

In this paper we have studied a two-species model with open boundaries and 
the conditions under which an invariant multiple shocks measure exists for it. 
We have found that from the MPF point of view the quadratic algebra of the 
model can be mapped to that obtained for the PASEP, if one defines the total 
density of particles, and also redefines the boundary conditions (the sum of 
the densities of particles at each site is defined as the total density of parti- 
cles at that site). Then we have been able to introduce a finite-dimensional 
Fock representation for this quadratic algebra, provided that the microscopic 
rates satisfy some constraints. It terms of the total density of particles these 
constraints are nothing but the necessary conditions for the existence of mul- 
tiple shocks in the system. The associated quadratic algebra of the model has 
also a one-dimensional representation under some restrictions which are the 
conditions for the existence of a uniform stationary product measure for the 
model. For the single-shock case our results obtained from MPF analysis agree 
with those obtained in [17] from study of the product shock measure dynamics. 
In [12] the authors have shown that for three different families of one-dimensional 
driven diffusive models with open boundaries, the stationary state distribu- 
tion function can be written as a linear superposition of single shock measures, 
provided that some constraints are fulfilled. In [18] the author has studied the 
same models and shown that from the MPF point of view these constraints 
are the necessary conditions for the existence of two-dimensional representa- 
tions for their quadratic algebras. The problem now worth studying is whether 
every finite-dimensional representation of the quadratic algebra of a given one- 
dimensional driven diffusive model with open boundaries is associated with 
expression of its stationary state distribution function in terms of superposi- 
tion of product shock measures. This is still under investigation. 
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